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Abstract 

In a homogeneous pionic gas system, chemical nonequilibrium pro- 
cess is considered to understand its effect in the expansion processes 
that are realized immediately after heavy ion collisions. The chemical 
relaxation time is calculated by incorporating the tt+tt «-> Tr+ir+ir+ir 
reaction, which is given in the second order of perturbation in the a 
model. The 7r + 7r <-» 7r + 7r + 7r + 7r reaction is assumed to be less 
frequent than the tt + tt <-> tt + tt scattering that is expected to estab- 
lish the local equilibrium, and hydrodynamical equation is solved for 
various initial conditions. It is shown that the relaxation time is of the 
order of lOOfm and does not have a significant effect on the expansion 
process, which implies that the pion number freezeout takes place at 
an early stage of the expansion. 
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1 Introduction 



It is widely believed that Quantum Chromo Dynamics (QCD) is 
the underlying theory of strong interaction and that QCD exhibits 
confinement-deconfinement phase transition at high temperature (~200 
MeV[[l|].) Since understanding the property of this phase transition is 
crucial for hadron physics, many efforts have been made both exper- 
imentally and theoretically [Q]. In particular, heavy ion reaction is 
expected to be a useful method to establish high temperature with 
free quarks and gluons and to investigate the process of hadroniza- 
tion. 

Heavy ion collision is considered to undergo several steps shown in 
Fig.L: 

1. Two highly Lorentz contracted nuclei collide with each other. 

2. They pass through each other and a hot system is formed in 
between, where quarks and gluons are excited @. 

3. Frequent collisions and reactions among the excited elements 
lead to establishing local equilibrium. Therefore this state can 
be treated by hydrodynamics Q ||. 

4. As the system expands, cooling and subsequent phase transition 
process proceeds Q. 

5. A lot of pions ( and other light hadrons) are created under the 
process of the phase transition ||. 

6. When the phase transition terminates, there exists hadronic gas 
which can also be treated by hydrodynamics. Since the system is 
hot, the collisions and reactions between the hadrons are frequent 
and the total number of hadrons is not conserved. 

7. As the cooling goes on, the reaction to change the total number 
of hadrons becomes unlikely. At this stage, the multiplicity of 
pion (and other light hadrons) is fixed (freezeout of number) 

8. When the expansion proceeds further and the mean distance be- 
tween pions is equivalent to the average variation of macroscopic 
quantity such as temperature, energy density and pressure, the 
hydrodynamical picture is not a good one any longer. Instead, 
the system can be considered to be composed of free pions with- 



2 



out any interaction. This stage is called the freezeout of temper- 
ature. 



Fig.l 

Therefore, the multiplicity of pion is determined at stage [?] not 
at the time when the phase transition terminates. Also, the energy 
distribution of pions is given by the thermal distribution at stage § 
In other words, hadronic final state of heavy ion collision does not 
directly reflect the property of the confinement-deconfinement phase 
transition 

Due to this feature, it has been proposed to use electromagnetic 
and weak probe II, which is free from the hadronic final state inter- 
action, to investigate the hadronic property of the phase transition. 
However, if we can understand the process of the final state interac- 
tion, it is possible to extract the information on the phase transition 
by hadronic probes. In particular, pion multiplicity can be used to 
obtain the pion number and the entropy immediately after the phase 
transition: if one knows the elementary chemical process to change 
the pion multiplicity such as tt + tt «-> tt + tt + tt + tt, one can see 
how the the pion multiplicity varies through the cooling process, and 
from the information on the pion multiplicity obtained in heavy ion 
reactions one can track the process backward up to the time when the 
phase transition is over. 

Whether the chemical reaction largely affects the cooling process 
or not depends on the balance between the cooling by the expansion 
and the speed of the chemical reaction. If the typical time scale of 
the chemical reaction is longer enough than that of the expansion, the 
freezeout of pion multiplicity takes place at an early stage of the ex- 
pansion and one can identify the final pion multiplicity with the pion 
multiplicity immediately after the phase transition in the first approx- 
imation. In this case we can estimate the entropy density immediately 
after the phase transition (s) by 

s = 4n, n : observed pion multiplicity, (1) 

which is obtained for massless classical particle systems with the par- 
ticle numbers conserved Q. If the former is shorter than the latter, 
the chemical reaction remains dominant through the expansion. In 
this case, the final pion multiplicity is quite different from the initial 
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one at the end of the phase transition, and Eq.(^) can not be used to 
obtain the entropy density. Instead, one must track the irreversible 
process backward, estimate the pion multiplicity at the end of the 
phase transition and then use Eq.(|l]). 

In this paper, we aim to obtain the speed of the chemical reaction 
7r + 7r *-* 7r + 7r + 7r + 7r based on the a model. In this model, the inter- 
action in the leading order of perturbation just gives it + tt *-* tt + tt 
scattering, which is expected to establish local equilibrium in a pionic 
gas. In the second order of the perturbation, ir + ir «-> tt + tt + tt + tt re- 
action is derived. If the latter reaction is less frequent than the former 
one, we can resort to hydrodynamical equation to track the change of 
pion number multiplicity. In section ||, the evaluation of the frequency 
of these two reaction is made and the criterion to establish local equi- 
librium is given. Then, we present hydrodynamical equation in the 
classical approximation. In section [| the hydrodynamical equation is 
solved numerically in a homogeneous pion gas for various initial tem- 
peratures and chemical potentials. The typical time scale of chemical 
reaction (relaxation time) is obtained and compared to the expansion 
rate. The last section || is devoted to summary and discussing open 
problems. 

2 Elementary chemical reaction in the 
cr-model and the hydrodynamical equa- 
tion 

We take cr-model to describe the tt + tt^tt + tt + tt + tt reaction. 
The interaction part of the Lagrangian is given by 

£int = -^ 2 + ^ 2 -/. 2 ) 2 , (2) 

where tt = tt 1 ^ + 7r 2 T2 + tt 3 t^ and f n is the pion decay constant. After 
shifting the a field as a —* a+f n , one obtains — X/Al(iT 2 +2f n a+a 2 ) 2 as 
the potential term. This term gives tt + tt *-* tt + tt scattering depicted 
in Fig. 2(a). Summing all the diagrams in Fig. 2(a) at threshold pi = 
(m, 0,0,0), (i = 1,2,3,4) gives 

V m a 2 { 3 ) m CT 2 V m 2 -{2m) 2 ' K ) 
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where m a (to) is the a (tt) meson mass. Since m a 2 = A/ vr 2 /3, we get 

- - A 1 

3 + 3 1 _ (2m)2 ( ) 

as an effective 4tt coupling constant. For m a ps 600 MeV, m/m a ps 1/4 
and the effective coupling constant is A/9. For simplicity, we hereafter 
replace the sum of Fig. 1(a) by just one direct Air coupling with its 
coupling constant A/9. 

The tt + tt^tt + tt + tt + tt reaction is given in the second order 
of perturbation (Fig. 2(b)). In the second order, another diagram is 
possible for the tt+tt <-> tt+tt+tt+tt process (Fig.2(c)). This diagram, 
however, is not taken into account because the denominator of the 
propagator p 2 — m 2 is of the order of (3m) 2 — m 2 and the contribution 
of Fig.2(c) is 0(1/10) of Fig.2(b). All the possible diagrams of type (b) 
are shown in Fig.2(d). For simplicity, we replace Fig2.(d) by Fig2.(e), 
which means the effect of interference is not exactly taken into account. 

Fig-2 

In the next section, the hydro dynamical equation is solved as- 
suming local equilibrium that is expected to be established by the 
tt + tt <-> tt + tt scattering. Since tt + tt <-> tt + tt + tt + tt reaction 
drives the system away from equilibrium, tt + tt <-> tt + tt scattering 
must be frequent enough so that the disturbed nonequilibrium state 
can swiftly return to equilibrium state. 

Once local equilibrium is established, local temperature T(x) as 
well as local chemical potential fj,(x) is defined. The tt + tt <-> tt + tt 
actually means 7r°+7r° <-> tt + +tt~ and 7r i +7r J ' <-> 7r l +7r ? = +, — , 0). 
The first one gives the equilibrium condition 2/i°(x) = ^ + {x) + [i~(x). 
Since the charge conservation means l~i + {x) = jjT(x) in neutral pionic 
gas, we get /U + (x) = [i~{x) = fi°(x) or n 7r+ (x) = n 71 " (x) = n 71 " (x) 
(n: pion density.) This means that all the equations are identical for 

TT°,TT + and 7T~. 

In other words, all the calculation is made in the same way as in 
A0 4 theory except for a numerical factor stemming from the isospin 
degrees of freedom, because the essential part of the a model is only 
the (A/4!)vr 4 term. 

As we have mentioned before, the tt + tt <-> tt + tt scattering should 
be more frequent than tt + tt ^tt + tt + tt + tt reaction in order to have 
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the system locally equilibrated. The frequency of7r + 7r<->7r + 7r and 
7r + 7r^->7r + 7r + 7r + 7r process are calculated in Appendix A: 



I 



■K-\-n^n-\-n 



dpi dp 2 dp 5 dp4 

(2vr) 3 2pi° (2vr) 3 2p 2 (2vr) 3 2p 3 (2vr) 3 2p 4 

x^) 2 (2tt) 4 ^i +p 2 -ps - p 4 )e-^ + ^° +2 ^ 

A s n 1 



2 m 



d|p||p| V 1 "^ 3 ^) (5) 



7 = S 3 ( ^ )4 T28^F e " 4/3m4(e2/3 " " ^ r AX + ^{x)e^d x 

(6) 

where (3 = 1/T, T is the temperature, \i the chemical potential, m 
the pion mass, Pq the total energy of the colliding pions, p their total 
spatial momentum, Po/m = 4 + x and g{x) a function independent 
of temperature and chemical potential whose definition is given in 
Appendix A. The factor ||| comes from the isospin degrees of freedom 
(see Appendix B.) Since no further analytical calculation is possible, 
we carried out the integrals numerically for mj3 ~ 1 and fi/m ~ 1 to 
get a rough estimate for / 7r + 7r ^ 7r + 7r anc [ j, it was found that A -C 1000 
is required for / 7r + 7r ^ 7r + 7r ^> / Since A ~ 100 is taken for the ordinary 
ex-model, local equilibrium is established in the present case.[] 

Under local equilibrium, relativistic hydro dynamical equation [p~C|] 
is given as : 

dpN*" = 21 (7) 
d^ u = (8) 

for each species of pion, where N % ^{i = +, — , 0) is pion number current 
and T l ^ u is the energy momentum tensor for 7r + ,7r~ and 7r° respec- 
tively. As is mentioned above, these fundamental equations take the 

1 Since 7r + 7r^-i'7r + 7ris 0(X 2 ) and 7r + 7r^->7r + 7r + 7r + 7ris 0(A 4 ), it may seem 
that A -C 1 is required for local equilibrium. However, numerical factors coming from the 
phase space integral such as l/(2n) 3 for each external line and e _4/3 in Eq.(||) suppress 
the frequency of the 7r + 7r<^7r + 7r + 7r + 7r process. Thus A ~ 100 satisfies the criterion 
for local equilibrium even if A is bigger than 1. 
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same form for all the species of pions. Hence, we hereafter suppress 
the indices % = +, — ,0. Eq.(|7|) describes the change of the total pion 
number and is called rate equation, whose right hand side 21 gives 
the rate of particle production and annihilation. The factor 2 in the 
right hand side comes from the fact that two pions are created for 
one elementary reaction in Fig. 2(b). Eq.(||) is the energy-momentum 
conservation identity. 

and T^ v are expressed in terms of the pion distribution function 

/: 

J (2vr)V W 



J (2vr)3 p o J - {W) 

3 Numerical solution for the hydro- 
dynamical equation and chemical relax- 
ation time in a homogeneous pionic gas 

In this section, we consider homogeneous pionic gases to obtain 
the typical time scale of the chemical reaction tt + n ^ it + n. In a ho- 
mogeneous gas, the hydrodynamical equation Eqs.([?])(||) is drastically 
simplified toQ 



§-»■ « 12 » 

where n = N°(e = T 00 ) is the pion density (energy density). To 
express and T^ v in terms of the local temperature T and local 
chemical potential n, we need to have the explicit form of the distri- 
bution function /. In this paper, we take classical distribution 

f(p, x, t) = exp{-/3(x, t)(^/p 2 + m 2 - fx(x, t))}. (13) 
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Since we investigate homogeneous pionic gas, there is no x dependence 
in any macroscopic quantity and T and \x are the functions of time t 
alone. Then, we get 

nit) = ^ /T Tm*K 2 Q (14) 



e(t) = ^V{3T*K 2 (^) + TmK^)} (15) 



Combining Eqs.(|ll|)(|l2|) with Eqs.(14)(15), one gets a set of differen- 
tial equation for T(t) and n(t). The energy conservation law Eq.(|l2|) 
means e(t) = e(t = 0) giving the equation to determine [i{t) from T{t) 
and initial conditions To = T(0),/j,q = /x(0): 

„ = —T + T \n{{^f K2{m,T f' \ + }, (16) 

^ T XK T J 3K 2 {m/T) +m/TK l {m/T) V ' 

with the abbreviation /i = /x(i),T = T(t). Also from de/dt = 0, we 
can express d\x/dt in terms of T{t), n{t),dT / dt,T$ and fiQ,: 

A = (| - A)T (17) 



12K 2 (m/r) + bjm/T^im/T) + (m 2 /r 2 )Er (m/r) 

3K 2 (m/T) + (m/T)Kx(m/T) ' 1 j 

where /i (T) means dfi/dt (dT/dt). By eliminating ^z, and (i and re- 
defining all the quantities to be dimensionless, i.e. T/m — > T = 
1/(3, fi/m — > /i, A 4 mi — ► i, we can rewrite the hydrodynamical equa- 
tion as 

{(3-^)K 2 (/3)+/3K 1 (/?)}T = 2//A 4 . (19) 

Since I is proportional to A 4 (see Eq.@), the A 4 term in the redefini- 
tion of t makes the equation independent of A. 

By solving this differential equation, one gets the t dependence of 
T and from Eq.(|ilf) fi(t) is also obtained. The results of numerical 
calculation are shown in Figs. 3, 4, 5, 6 for initial conditions (To,/^o) = 
(0.5, -1), (1.5, -1), (0.5, 1), (1.5, 1) respectively. 

Figs. 3,4,5,6 
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From these figures, one can see that the relaxation is faster for 
hotter and denser (larger /^o) systems, which is reasonable consider- 
ing that the chemical reaction is more frequent for hotter and denser 
systems. Since the relaxation is not exactly exponential, we define re- 
laxation time by ^(t re iax) = W)/ e - The relaxation times with A = 120 
are shown in Fig.7 for various initial conditions —1 < /io < 1,0.5 < 
To < 1.5. For a dense and hot system (e.g. /Uo = 1 and To = 1.5), 
^reiax is of the order of lOOfm. This time scale is expected to be larger 
than the typical time scale of expansion. Thus, we have found that 
the 7r + 7r^7r + 7r + 7r + 7r reaction does not play an important role 
in an expanding pionic gas system. In other words, the pion number 
freezes out at an early stage of the final state of a heavy ion reaction. 
This implies that the hadronic probe is a good one to extract informa- 
tion on the phase transition and that s = 4n is an appropriate one to 
evaluate the entropy density immediately after the phase transition, 

Fig.7 



4 Summary and open problems 

We have investigated the effect of7r + 7r <-> tt + it + tt + tt chemical 
reaction in the framework of hydrodynamics. To understand the typ- 
ical time scale of the tt + tt^tt + tt + tt + tt reaction, we calculated 
the relaxation times in homogeneous pionic gas systems. It turned 
out that the relaxation time is of the order of lOOfm at shortest and 
is longer than typical time scale of the expansion of the pionic gas. 

Although we found that the chemical tt+tt ^ tt+tt+tt+tt reaction 
does not play an important role within the framework of the <r-model, 
there still are a lot of things to be done. The first one is to take into 
account the effect of the baryon in the central region, which has been 
neglected so far. In the real case, thermally excited NN pairs exist 
in the central region. Since the ttN interaction is much larger than 
the tttt interaction, the processes such as./V-|-7r<->./V-|-7r-|-7r might 
be significant even if the number of the excited NN pairs is small. 
Furthermore, if there is a stopping in heavy ion collisions and there 
exists baryon number coming from the colliding nuclei, there would 
be a big change of pion multiplicity. 

Secondly, to make reliable quantitative discussion on the role of the 
tt + tt^tt + tt + tt + tt reaction, the calculation must be made without 
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classical approximation incorporating the effect of interference. This 
is expected to make the collision less frequent and lead to a larger 
relaxation time, keeping the essential statement in the former section. 
The third thing is to take the medium effect into account. So far, 
we have fixed the parameters in the a-model: the coupling constant 
A and the pion mass m. It is known, however, that these quantities 
depend on temperature and chemical potential. The correction due 
to this effect might not be negligible. 

Even if some of these effects lead to a smaller relaxation time 
and make the effect of the change of the pion number a little more 
important, we still do not expect them to change our results qualita- 
tively because of the following reason. In an expanding system, the 
chemical potential at the end of the phase transition is taken to be 
0, because there are strong interaction between QGP and pions and 
pion can change its multiplicity easily. Therefore, at an early stage 
of the expansion in which the system is hot enough, the chemical po- 
tential is close to and the system is close to chemical equilibrium. 
So, ir + ir^ir + ir + TT + Tr reaction etc. can not play an important 
role. As the system expands, the chemical potential deviates from 0. 
However, the system is not hot enough and the relaxation time is too 
long to make the effect of"7r + 7r^7r-|-7r-|-7r-|-7r reaction significant. 
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Appendix A 



In this appendix, we show how to calculate the frequency of it + ir <-> 
7r + 7T scattering and 7r + 7r^7r + 7r + 7r + 7r reaction. For simplicity, 
we neglect the isospin degrees of freedom, which will be taken into 
account in Appendix B. 

Under classical approximation, the frequency of the process de- 
picted in Fig. 2 (a) per unit volume and unit time is given by 



dpi dp2 dps dpi 

(27r)32 Pl o (27r)32p 2 ° (2vr)32p 4 



x (^) 2 (27r) 4 «5( Pl + p 2 -p 3 - P4 )e-^ 1+ ^° +2 ^A.l) 

We first carry out the integral over dps and dp4 . Rewriting the integral 
in an explicitly Lorentz invariant form by using 

/ (2^%? = / dP3 6{P ' ] 5{PI ~ ^ (A - 2) 
and then changing the variables to p = ps + P4, q = we get 



f dp 3 dp 4 x , . , 

J *(P3+P4-P1-P2) 



Pl /i 4m2 



= y dp dp J ^-^dqo j^2tt 2 J d \q_\\q\ 2 -8(p- Pl -p 2 ) 
x y d(cosfl)5( 



a P°1 \Xf\ I / 2 , P 2 2\ 

cos9 ~ rinWw -\% + -r- m )- 



p q V 4 



The integral over dp\ and dp 2 is made essentially in the same way. 
However, the factor J\ — , 4 ™ 2 ^ in the above equation now prevents 



O1+P2) 2 

full analytical calculation. The final result is 



V 32 (2tt) 5 i 2m 

V 7 ^ 2 " 4 ™ 2 2 £ 4^2 



'" p||p| '"-(n^y (A - 4) 
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where Pq = (p% + p 2 )° and p = Pi + P2- 

The collision frequency for Fig. 2(b) is made almost in the same 
manner. All the possible diagram of type Fig2.(b) is given in Fig2.(d). 
We approximate Fig2(d) by Fig2.(e), which means that we have re- 
placed the interference term by Fig. 2(e). This approximation makes 
it possible to carry out the integral over the relative variables between 
pi and P2 and between p 3 and p 4 . The rate of the particle production 
of the processes 7r + 7r< — 7r + 7r + vr + 7r depicted in Fig.2(e)is written 
as 

1 f dpi dp 2 dp 3 dp4 dp 5 dp 6 

4!2! J (2^) 3 2p!° (2vr) 3 2p 2 ° (2vr) 3 2p 3 (2vr) 3 2p 4 (27r) 3 2p 5 ° (2^) 3 2p 6 ° 
x e 2l3 > 1 e- f3{p5+P(i)0 \6M\ 2 (2iT) 4 5(p 1 +p 2 + P3 + P4 - P5 -p (i ), (A.5) 

where M is the matrix element of the Fig. 2(b): 

M = d? 2 , \ - 2 - (A.6) 

9 m z - (pi +P2 - P 5) 

The rate for the reverse reaction ir + ir + ir + n <— ir + ir is given by re- 
placing e^e-^+P 6 ) by e 4 /?M e -/3(pi+P2+P3+P4) . The difference from 



Eq.( A.l ) is that the propagator in M hampers analytically carrying 
out the integral of d(p\ + p 2 ) and d(p 3 + p 4 ) as is shown below. 

Taking P1+P2 = p, Pl ~ P2 = q,p 3 + p 4 = p', P3 ~ P4 = q' and carrying 



out the integral over dq and dq' as is done in Eq.( [A2| ), we get 
1 




/ 4m z f 

, dp\ l -j- / dp 5 dp 6 

P z Jp> 2 >4m? V P ■> 

xS( P5 2 - m 2 )5( P6 2 - m 2 )9(p 5 )6(p 6 °)5(p+p f + P5 + P(i ) 

m ~ \P ~ P5) 

Then, we rewrite ps + p$ = pt, P5 ~ P6 = q t: p + p' = P, 2^£L = Q and 
define 

P = (P°,0,0,|P|) 

Q = (Qo) IQI sin^i, 0, |Q| cos^i) 

It = {qt°, |q| sin 9 2 cos <p, |q| sin Q 2 sine/?, |q| cos 62). (A. 8) 
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Pt can be eliminated due to the total momentum conservation 5 func- 
tion. The other two 5 functions are written as 



5(p 5 2 - m 2 ) 5(p 5 2 - m 2 ) = 8(— + Q 2 - m 2 ) S(2P ■ Q) (A.9) 

and enable us to eliminate cos 62 and |qd. The integration over dip 
is easily made and after some manipulation to make the expression 
convenient for Monte Carlo calculation, we get 



/ = 



,A 



128(2vr) 9V 9 



/ x 2 (x + 8)g{x)e- (3x dx, 
Jo 



(A.10) 



rl rl rl rl rl 

g(x) = / dx 2 / (ix 3 / dx4 / dz dy 

Jo J-i Jo J-i J-i 



x 6{— + Q 2 -A- \P-Q\) x 2 x\ 
with dimensionless variables redefined by 



(A.11) 



D Pl + P2 + P3 + P4 P5+P6 ^ (pi +P2-P3- Pa) (P5 ~ P&) 

r — — M — - , qt - 



m 



m 



2m 



2m 



P° = 4 + x, |P| = V P° 2 ~ 16x 2 = yjx{x + 8)x 2 , Qo = -x 3 , |Q| = lx 4 , 



Qt 
I = 



lp 2 - |P| 2 



1 + (qt ) 2 , 



I pO 2 _ 1 R p2 

+ Q° , a = — -Q 2 + 2Q°q t ° - 2|Q||q t |;^ 2 , 



P°<k\ 



b = 2|Q||q t |^/l — z\\j\ — z 2 , z = z\ = cos6±,Z2 = cos9 2 = j—fp— |- (A. 12) 



|p||q* 



The 9 function comes from the constraints 



p 2 = (_ + Q) 2 > Am 2 and p' 2 = ( - - Q) 2 > Am 2 (A.13) 
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The advantage of taking the classical distribution is that the thermal 
factor is given as e _/3P ° and all the variable other than the total energy 
can be integrated regardless of temperature. The temperature inde- 
pendent part is written as g(x) whose x (or Pq) dependence is shown 
in Fig.A. 

Fig.A 



Appendix B 

In this appendix, we consider the effect of the isospin degrees 
of freedom. For simplicity, we take tt ,^ 1 ,^ 2 not ir°,ir + ,ir~ as three 
species of pions. One can easily confirm that the calculation based on 
7r°,7r + ,7r _ leads to the same conclusion. 

The factor attached to four point vertices depends on whether the 
coupling concerns only one species of pions or not: —A for the 4tt° 
(47T 1 and 4vr 2 ) vertices and -A/3 for the 27r -2vr 1 (2vr 1 -27r 2 and 2vr 2 - 
27T°) vertices (see Fig.B.l.) 

Fig.B.l 

First, let us take the ir° + ir° <-> 4tt processes with every possible 
combination of 47r such as 7r + +7r — + 7r° +7r°. All the possible Feynman 
diagrams are shown in Fig.B.2. Since the three (two) diagrams in 
Fig. B. 2(c) and Fig.B.2(c') (Fig. B. 2(d)) have the same initial and final 
state, they must be added before taking the square and interfere with 
each other. Then we have to take possible combinations of allocating 
the external lines to the initial and final pions. As is mentioned in 
Chap. 2, we make the approximation that all the kinds of phase space 
integrals give the same result. So, we just multiply the number of 
combinations to the simplest diagrams shown in Fig.B.2 and then put 
appropriate factors to avoid overcounting of the initial and final states. 
Now that phase space integrals are the same for all the diagrams, we 
just need to calculate the factors coming from the vertices (—A and 
—A/3) and from the combination of allocating the external lines shown 
above. 

Fig.B.2 

These calculations are easily made and the diagrams in Fig.B.2 
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give |A 4 + |A 4 for tt 1 + tt 1 <-> 47r. The calculations for 7T 1 + ir 1 <-> 4-7T 
and 7T 2 + 7r 2 4-7T are made exactly in the same way, and one obtains 
3(|A 4 + jrA 4 ) as the sum of tt 1 + tt 1 <-> 47r (i = 1, 2, 3) reactions. 

Next, we need to consider those reactions from (to) two different 
species of pions to (from) all the possible four pions. As an example, 
we take ir° + tt 1 <-> Air. The possible processes are depicted in Fig.B.3. 
The calculations are essentially the same as before, and we get |y A 4 
as the sum of these diagrams. The same results are obtained for the 
other two processes tt 1 + tt 2 «-> 4tt and tt 2 + tt° <-> 4-7T. 

Fig.B.3 

Thus all the contribution of isospin degrees of freedom is 3(| + 
| + |f )A 4 , or the increase and decrease for each species of pion is 
characterized by (| + | + ff )A 4 . In other word, the effect of isospin 
degrees of freedom can be taken into account by replacing |A 4 , which 
is the factor in the case of 7r° + tt° <-> ir° + 7r° + 7r° + 7r° given by 
Fig.B.2(a), by (f + § + f )A 4 . The ratio (| + i + f )A 4 /f A 4 is the 
numerical factor ||| in Eq.(^). 
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Figure captions 

Fig.l: The schematic picture of heavy ion collision in t-z space. 

Fig. 2: The elementary process under consideration. (a):7r+7r <-> tt+tt 
scattering. (b):7r + 7r <-> 7r + 7r + 7r + 7r process. (c):the other pos- 
sible 7T + 7r^7r + 7r + 7r + 7r process in the second order of per- 
turbation in the a model. This diagram is neglected. (d):All the 
possible detailed diagram of type (b). (e): Approximated diagram 
for (d) used in the numerical calculation. All the 4-zr vertices in 
(b),(c),(d),(e) are the effective coupling obtained from the sum 
of diagrams in (a). 

Fig.3.4.5.6: The time dependence of (a) temperature T, (b):chemical 
potential /j,, and (c)dog \fx\ for initial conditions (To, no) = (0.5, —1.0) 
(Fig.3), (1.5,-1) (Fig.4), (0.5,1) (Fig.5) and (1.5,1.0) (Fig.6). The 
temperature T, chemical potential fi, time t are all dimensionless 
ones defined in the text. 

Fig.7: log 10 trelax for various initial conditions with A = 120. 

Fig. A.: x dependence of the temperature independent function g(x). 

Fig.B.l: Two kinds of 47r coupling in the a model. 

Fig.B.2: All the possible Feynman diagrams for 7r° + tt° <-> Air 

Fig.B.2: All the possible Feynman diagrams for ir° + tt 1 <-> 4-zr 
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